Dynamics and constraints of the Unified Dark Matter fiat cosmologies 



Spyros Basilakos 1 and Georgios Lukes-Gerakopoulos 1 ' 2 

1 Academy of Athens, Research Center for Astronomy and Applied Mathematics, 
Soranou Efesiou 4, GR- 11527, Athens, GREECE 
2 University of Athens, Department of Physics, Section of Astrophysics, Astronomy and Mechanics 

We study the dynamics of the scalar field FLRW flat cosmological models within the framework 
of the Unified Dark Matter (UDM) scenario. In this model we find that the main cosmological 
functions such as the scale factor of the Universe, the scalar field, the Hubble flow and the equa- 
tion of state parameter are defined in terms of hyperbolic functions. These analytical solutions 
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large and small scale dynamics which are in fair agreement with the predictions by the above dark 
energy models although there are some differences especially at high redshifts. 
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1. INTRODUCTION 

The detailed analysis of the available high quality cos- 
mological data (Type la supernovae [l[ , [1] ; CMB [l| , 
[H, etc.) leads to the conclusion that we live in a flat 
and accelerating universe. In order to investigate the 
cosmic history of the observed universe, we have to in- 
troduce a general cosmological model which contains cold 
dark matter to explain the large scale structure cluster- 
ing and an extra component with negative pressure, the 
vacuum energy (or in a more general setting the "dark 
energy"), to explain the observed accelerated cosmic ex- 
pansion (Refs. [H, i, i, H| and references therein). The 
nature of the dark energy is one of the most fundamental 
and difficult problems in physics and cosmology. There 
are many theoretical speculations regarding the physics 
of the above exotic dark energy, such as a cosmologi- 
cal constant (vacuum), quintessence, k— essence, vector 
fields, phanton^ tachyons, Chaplygin gas and the list 
goes on (see 0, |, 0, [1,1, Hfl HO, M M, M EE M M M 
and references therein). 

Such studies are based on the general assumption 
that the real scalar field <fi rolls down the potential 
V(6) and therefore it could resemble the dark energy 
[1, S El EI US HI, El . This is very important because 
scalar fields could provide possible solutions to the cos- 
mological coincidence problem. In this framework, the 
corresponding stress-energy tensor takes the form of a 
perfect fluid, with density p^ = <p 2 /2 + V{4>) and pres- 
sure Ptf, — 2 /2 — V ((/)). From a cosmological point of 
view, if the scalar field varies slowly with time, so that 
(/> 2 /2V <C 1, then w = P^/p^, ~ —1, which means that 
the scalar field evolves like a vacuum energy. Of course 
in order to investigate the overall dynamics we need to 
define the functional form of the potential energy. The 
simplest example found in the literature is a scalar field 



with V{4>) oc 4> 2 (see for review [Tl|, (23[) and it has been 
shown that the time evolution of this scalar field is dom- 
inated by oscillations around = 0. Of course, the issue 
of the potential energy has a long history in scalar field 
cosmology (see [24|, |25j, |26|, |27|, |28[ and references therein) 
and indeed several parameterizations have been proposed 
(exponential, power law, hyperbolic etc). 

The aim of the present work is to investigate the obser- 
vational consequences of the overall dynamics of a family 
of flat cosmological models by using a hyperbolic scalar 
field potential which appears to act both as dark matter 
and dark energy [2t|. To do so, we use the traditional 
Hamiltonian approach. In fact, the idea to build cos- 
mological models in which the dark energy component 
is somehow linked with the dark matter is not new in 
this kind of studies. Recently, alternative approaches to 
the unification of dark energy and dark matter have been 
proposed in the framework of the generalized Chaplygin 
gas [30|, HH and in the context of supersymmetry . 

The structure of the paper is as follows. The basic the- 
oretical elements of the problem are presented in section 
2 by solving analytically [for spatially flat Unified Dark 
Matter (UDM) scalar field models] the equations of mo- 
tion. In section 3, we present the functional forms of the 
basic cosmological functions [a(t), (p(t) and H(t)]. In sec- 
tion 4 we place constraints on the main parameters of our 
model by performing a joint likelihood analysis utilizing 
the SNIa data 0] and the observed Baryonic Acoustic Os- 
cillations (BAO) [HI and [H[ ■ In particular, we find that 
the matter density at the present time is tt m ~ 0.25 while 
the corresponding scalar field is <fio ~ 0.42 in geometrical 
units (0.084 in Planck units). Section 5 outlines the evo- 
lution of matter perturbations in the UDM model. Also 
we compare the theoretical predictions provided by the 
UDM scenario with those found by three different type 
of dark energy models namely A cosmology, parametric 
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dark energy model and variable Chaplygin gas. We verify 
that at late times (after the inflection point) the dynam- 
ics of the UDM scalar model is in a good agreement, 
with those predicted by the above dark energy models 
although there are some differences especially at early 
epochs: (i) the UDM equation of state parameter takes 
positive values at large redshifts, (ii) it behaves well with 
respect to the cosmic coincidence problem, and (hi) be- 
fore the inflection point the cosmic expansion in the UDM 
model is much more decelerated than in the other three 
dark energy models implies that the large scale struc- 
tures (such as galaxy clusters) are more bound systems 
with respect to those cosmic structures which produced 
by the other three dark energy models. Finally, we draw 
our conclusions in section 6. 



ANALYTICAL SOLUTIONS IN THE FLAT 
SCALAR FIELD COSMOLOGY 



the scalar field potential energy V(4>), which is not an 
easy task to do. Indeed, in the literature, due to the un- 
known nature of the dark energy, there are many forms 
of potentials proposed by several authors (for a review 



see [T3|, 16 1, [38|) which describe differently the physics 
of the scalar field. It is worth pointing out that for some 
special cases analytical solutions have been found (Refs. 
[2^, [2^, [23, [H, [la [4l|| and references therein). As an 
example, if the potential V((j>) is modeled as a power law 
4> n , then the energy density of the scalar field evolves like 
p,f, cx a -6n/(n+2) means that, for n = 2 or n = 4 

the corresponding energy density behaves either like non 
relativistic or relativistic matter. In this work, we have 
used a functional form of V(<f>) (see [42j) for which we 
solve the previous dynamical problem analytically. This 
potential corresponds to the so called Unified Dark Mat- 
ter (hereafter UDM) scenario [H|, [IH, [|||: 

V((j)) = Cl COsh 2 (V <j)) + C 2 CjCi, C2 6 R , (5) 



Within the framework of homogeneous and isotropic 
scalar field cosmologies it can be proved (see Ref. [351 ]) 
that the main cosmological equations (the so called 
Fricdmann-Lemaitre equations) can be obtained by a La- 
grangian formulation: 



-3aa 



v{4>) 



Ska 



(1) 



where a{t) is the scale factor of the universe, 4>(t) is the 
scalar field, V(4>) is the potential energy and k(— —1, 0, 1) 
is the spatial curvature. Indeed the equations of motion 
take the following forms [70j: 



v(4>) 



and 



(t> + 3~<t> + V'(4>) =0 
a 



(2) 
(3) 

(4) 



where the over-dot denotes derivatives with respect to 
time while prime denotes derivatives with respect to <f>. 
We would like to stress here that in this work we consider 
a spatially homogeneous scalar field 0, ignoring the pos- 
sible coupling to other fields and quantum-mechanical 
effects. On the other hand, introducing in the global 
dynamics a new degree of freedom, in a form of the 
scalar field 0, it is possible to make the vacuum energy 
a function of time (see [2(|, [3(|, |13|)- Note of course 
that the geometry of the space-time is described by the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) line el- 
ement. 

In order to study the above system of differential equa- 
tions we need to define explicitly the functional form of 



Following the Bertacca et al. 29] nomenclature, the 
real constants in eq. §5$ are selected such as c\ — Ci > 0. 
As expected there is one minimum at the point <f> = 0, 
which reads 



V min = V(0) = c 1 +c 2 



(0) 



We would like to point out that as long as the scalar 
field is taking negative and large values the UDM model 
has the attractive feature due to V(4>) cx e -227 ^ (24[. 
This property simply says that the energy density in <p 
tracks [41j the radiation (matter) component. In fact 
the UDM potential was designed to mimic both the dark 
matter and the dark energy. Indeed, performing a Taylor 
expansion to the potential around its minimum we get, 



v{4>) = v min + Cl v 2 4> 2 



ciV 4 



(7) 



which means that at an early enough epoch the " cosmic" 
fluid behaves like radiation [44j (V((j>) cx </> 4 ), then evolves 
to the matter epoch (V(cf>) cx (f> 2 ) and finishes with a 
phase that looks like a cosmological constant (see also 

Changing now the variables from (a, <f>) to (x\, x 2 ) us- 
ing the relations: 



Xl 



3'2 




3/2 sinh (V 
a 3/2 cosh (V 



(8) 



with T> 2 =3/8 the Lagrangian {1} is written: 
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3 3 
{x\ 2 + -c 2 xj) - (x 2 2 + -(ci + c 2 ) xl) 



1 



3 4 / 3 k(x 2 - x 2 ) 1 ' 3 



(9) 
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The scale factor (a > 0) in the UDM scenario is now 
given by: 



3(x 



1/3 



(10) 



which means that the new variables have to satisfy the 
following inequality: X2 > \x\\. 

It is straightforward now from the Lagrangian ^ to 
write the corresponding Hamiltonian: 



H 



1 



[(P 2 X1 ~ u\x\) - (p 2 X2 - uj^xD] 



(11) 



where p Xl = x\, p X2 = 
menta and cuf = |c2, w| — j 
"frequencies" with units of inverse of time and 



x\ denote the canonical mo- 
3 (ci + c 2 ) are the oscillators' 



a; 3 



C2 



(12) 




The dynamics of the closed FRLW scalar field cosmolo- 
gies has been investigated thoroughly in [42]. In partic- 
ular, for a semi-flat geometry (k — > 0) we have revealed 
cases where the dynamics of the system (see section 3.1 
in [42], orbit 5 in Fig. 1 scale factor vs. time and Fig. 4) 
is close to the concordance A-cosmology, despite the fact 
that for the semi flat UDM model there is a strong in- 
dication of a chaotic behavior. In this paper we would 
like to investigate the potential of a spatially flat UDM 
scenario (k = 0) since the analysis of the Cosmic Mi- 
crowave Background (CMB) anisotropics have strongly 
suggested that the spatial geometry of the universe is flat 
Technically speaking, in the new coordinate system 
our dynamical problem is described well by two indepen- 
dent hyperbolic oscillators and thus the system is fully 
integrable. Indeed in the new coordinate system the cor- 
responding equations of motion can be written as 



dH 

dx\ 



ijj\x\ 



Px 2 



dH 

8X2 



-<J\xi 



and it is routine to perform the integration to find the 
analytical solutions: 



x\ it) = A\ sinh {u>\t + 9\) 
x%{t) — A2 sinh (u2t + 62) 



(13) 



where A%, A 2l 9± and 62 are the integration constants 
of the problem. With the aid of eq. (fT3|) and assuming 
that the total energy (Ti) of the system is zero the above 
constants satisfy the following restriction: 



A\J\ 



AIlj 



a\ ,4 

= — = — 



(14) 



As expected the phase space of the current dynamical 
problem is simply described by two hyperbolas — 



2J2 



LO'X 



ujfAf whose axes have a ratio I/oji (i = 1,2). 



FIG. 1: Likelihood contours in the (Q m , \@i\) plane. The con- 
tours correspond to la, 2a and 3a confidence levels. The thick 
(thin) contours correspond to the SNIa (BAOs) likelihoods 
while the solid square is the best fit solution: Q m — 0.25 and 
6*i ~ —0.39. Insert Panel: The solutions within la contours 
in the (Q, m ,4>o) plane, where (f>o is the present value of the 
scalar field in geometrical units (8irG = c = 1). Note that for 
Planck units ( G = c = 1) we have to multiply the <^o-axis 
with (871-) -1 / 2 . Using the best fit solution we find <f>o — 0.42 
or 0.084 in Planck units. 



3. THE EVOLUTION OF THE UDM 
COSMOLOGICAL FUNCTIONS 

In this section, with the aid of the basic hyperbolic 
functions, we analytically derive the predicted time de- 
pendence of the main cosmological functions in the UDM 
cosmological model. 



3.1. Scalar field - potential versus time 

If we combine eq. (|13p together with the initial param- 
eterization [see eq.®] we immediately obtain the follow- 
ing expressions: 

nismh(uiit + 61) 



— = tanh(£></>) 

X2 



sinh(w 2 £ + 62) 



*(*) (15) 



and after some algebra the evolution of the scalar field 
becomes: 



6{t) = —In 



1 + tf (t) 



Using eqs. ([5]) and (fTrJ)) one can prove that 



1 - * 2 (t) 



(16) 



(17) 
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FIG. 2: Upper Panel: Comparison of the scale factor pro- 
vided by the UDM model (solid line) with the traditional A 
cosmology (short dashed line), VCG (dot dashed) and CPL 
(long dashed) dark energy models. Bottom Panel: The evo- 
lution of the scalar field. In the insert panel we present the 
behavior of the potential normalized to unity at the present 



time. Note, that to 
the universe. 



Now the range of ^-values for which the UDM scalar field 
is well defined [due to eg .([!$])] is: * e (-1, 1). Evidently, 
when the system reaches at the critical point = then 
^(t m ) — (or V m in = 4cj 2 k/3). For this to be the case 
we must have t m = —Q\ju\ and therefore, 9\ < 0. 



3.2. Scale factor - Hubble flow versus time 

Inserting eq. (fT"3")) into eq. (fTU|) the scale factor, normal- 
ized to unity at the present epoch, evolves in time as 



implies that the condition 

V{(f)i) -4>i 2 = 

should contain roots which are real and such that a S 
(0, 1). The above equation is solvable because for c\, c% > 
the potential energy [V(0)] takes only positive values. 
Knowing the integration constants (u>\, 9\, u>2,9 2 ), of the 
current dynamical problem we can calculate the inflection 
point by solving numerically the following equation: 

2* 2 (t,) -uI[k- <F 2 (t/)] [1 - <F 2 (t/)] = (20) 

with 

2 

y/nui 2J Vi smh(w 3 _^ + 63^) cosh(u>it + 0j) 

= — 2 (21) 

sinh (u 2 t + 9 2 ) 

where V\ = 1 and v 2 = —^/k. 

In addition, the Hubble function predicted by the 
UDM model can be viewed as the sum of basic hyper- 
bolic functions: 



o(t = - = 



sinh 2 (cj2< + 62) — Ksinh 2 {uj\t + 9\) 
sinh 2 (u> 2 to + 62) — Ksinh 2 (u>ito + 6\) 



1/3 

(18) 



where to is the present age of the universe in billion years. 
The constant 9\ is related to 9 2 because at the singular- 
ity (t = 0), the scale factor has to be exactly zero [see 
cq. lfTSll ]. After some algebra, we find that 



6» 2 = In ( -v^sinhfli + \J Ksinh 2 6q + lj > . (19) 

Furthermore, we investigate the circumstances under 
which an inflection point exists and therefore have an 
acceleration phase of the scale factor. This crucial period 
in the cosmic history corresponds to a(ti) = which 



_ a 1(x2X2 — x\X\) 2 



H(t) = - = ^ = ^ 2 f{t) = h ^{22) 



3(^2 X l) 



fit) 

7(*o) 



H 1 ~ 13.6Gyr is the present age of with 



fit) 



^ v% sinh {ujz-it + 63-i) cosh (u 3 _ii + 6» 3 _,) 

i=l 

sinh 2 {u>2t + 62) — K sinh 2 (wit + 81) 



-(23) 



where Hq is the Hubble constant. In this work we 
use H = lOO/iKms^Mpc" 1 with h = 0.72 [H| or 
Ho = /i/9.778 - 0.0736Gyr" 1 corresponding to t m 
H^ 1 ~ 13.6Gyr. Also we can relate the frequency u>2 
of the hyperbolic oscillator in the X2 axis with the well 
known cosmological parameters. Indeed, LO2 is given by 



0J2 = 



r. 



(24) 



while 9\ oc Hoto has no units. Notice that, VL m is the 
matter density at the present time. 



4. COSMOLOGICAL CONSTRAINTS AND 
PREDICTIONS 

In this work we use the so called Baryonic Acous- 
tic Oscillations (BAOs) in order to constrain the cur- 
rent cosmological models. BAOs are produced by pres- 
sure (acoustic) waves in the photon-baryon plasma in the 
early universe, generated by dark matter overdensities. 
First evidence of this excess was recently found in the 
clustering properties of the luminous SDSS red-galaxies 
(HI, [iH and it can provide a "standard ruler" with which 
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we can put constraints on the cosmological models. For 
a spatially flat FLRW we use the following estimator: 



[z 2 a H(a a )/H }W 



as y 2 H{y)/H a 



2/3 



(25) 



measured from the SDSS data to be A = 0.469 ± 0.017, 
where z s = 0.35 [or a 8 = (1 + z s ) _1 ~ 0.75]. Therefore, 
the corresponding Xbao function is simply written 



Xbao(p) = 



[A(p) - 0.469] 2 
0.017 2 



(26) 



where p is a vector containing the cosmological parame- 
ters that we want to fit. 

Also, we additionally utilize the sample of 192 super- 
novae of Davies et al. Q. In this case, the XgNia function 
becomes: 



1!)2 



XSNIi 



^fap) - H obB (Oi) 



(27) 




FIG. 3: Residual magnitudes (relative to the traditional 
Qa = 0.74, ttm = 0.26 model) of 192 SNIa data (open points) 
from 0] as a function of redshift. For comparison we plot the 
A(m — M)udm-a (solid line), A(m — M)cpl-a (long dashed 
line) and A(m — M)vcg-a (dot dashed line). 



(dof = 192) in a very good agreement with the 5 years 
WMAP data [3| . The inflection point takes place at 



where a% = (1 + z^) -1 is the observed scale factor of the 
Universe, Zi is the observed redshift, fi is the distance 
modulus \i = m — M = 51ogc?L + 25 and di,(a,p) is the 
luminosity distance 



d h (a,p) 



a J a V 2 H{y) 



(28) 



where c is the speed of light (= 1 here). Finally we 
can combine the above probes by using a joint likelihood 
analysis: 



Aot(p) = £bao x £ S Nia xL(p) = Xbao 



XsNIa 



in order to put even further constraints on the param- 
eter space used. Note, that we define the likelihood 
estimator [71j as: Cj oc exp[— x?/2]. 



4.1. The standard A-Cosmology 

Without wanting to appear too pedagogical, we remind 
the reader of some basic elements of the concordance A- 
cosmology. In this framework, the normalized scale factor 
of the universe is 



1/3 



sinh 2/3 u 2 t 



1-0 

The Hubble function is written as 

2 2 
H(t) = -CJ2 f{t) = -w 2 cothw 2 i 



(29) 



(30) 



Comparing the A model with the observational data (we 
sample tt m S [0.1, 1] in steps of 0.01) we find that the 
best fit value is Q m = 0.26 ± 0.01 with Xtoti^m) - 195 



1 



• sinh 









1/3 




4 = 


[2(i -n m )_ 





(31) 



Therefore, we estimate uo 2 ^ 1.29H ^ 0.095Gyr" 1 , tf ~ 
0.51io and af ~ 0.56. The deceleration parameter at the 
present time is qo = — a/d 2 | a= i ~ —0.61. 



4.2. The parametric Dark Energy model 

In this case we use a simple parameterization for the 
dark energy equation of state parameter which is based 
on a Taylor expansion around the present time (see 
Chevallier & Polarski [11] and Linder (43], hereafter CPL) 



w(a) = wo + wi(l — a) 
The Hubble parameter is given by: 



(32) 



H(a) = H a 



n, 



-3 



+ (1 - n m )a- 3(1+Wo+Wl) e awi( °- 1) 



1/2 



(33) 

where wo and wi are constants. We sample the un- 
known parameters as follows: wo G [—2, —0.4] and wi G 
[-2.6,2.6] in steps of 0.01. We find that for Q m = 0.26 
the overall likelihood function peaks at wo = — 1.20^ 20 
and wi = 1.141^'g while the corresponding x 2 ot (wo,wi) 
is 193.6 (dof = 191). The deceleration parameter at the 
present time is qo — —0.83. 



4.3. The Variable Chaplygin Gas as an alternative 
to dark energy 

Let us consider now a completely different model 
namely the variable Chaplygin gas (herafter VCG) which 



6 



corresponds to a Born-Infeld tachyon action [H, |49[ . Re- 
cently, an interesting family of Chaplygin gas models was 
found to be consistent with the current observational 
data [5(|. In the framework of a spatially flat FLRW 
metric, it can be shown that the Hubble function takes 
the following formula: 

1/2 



H(a) = H fl b a- 3 + (1 - Q b )y/B s a- 6 + (1 - B s )a~ 

(34) 

where Vt b ~ 0.021ft. 2 is the density parameter for the 
baryonic matter [5l| and B s G [0.01,0.51] in steps of 
0.01 and n G [—4,4] in steps of 0.02. The corresponding 
effective equation of state parameter w(a) = Pde/pde 
is related to H(a) according to 



w(a) 



1 - 



2 dlnH 

3 da 



1 



Ha\2 



(Ha 
\ H 



,-3 



(35) 



while the effective matter density parameter is: = 
fij, + (1 — Qb)-\/Ss- We find that the best fit parameters 
are B s = 0.07 ± 0.02 and n = 1.06 ± 0.33 (f2^ ~ 0.29) 



with Xtot{ B s 



193.7 (dof = 191) and the present 



value of the deceleration parameter is qo ~ —0.60. 



4.4. The UDM comparison with other Dark energy 
models 



In order to predict analytically the time evolution of 
the main cosmological functions [</>(i), a(t), H(t) and 
w(i)] we have to define the corresponding unknown con- 
stants of the problem (u>\, 0i,W2, 62)- At the same time, 
from the restrictions found in section 3 (see eqs. [12] fl9l 
and!24p. we can reduce the parameter space to (fi m ,#i). 
We do so by fitting the predictions of the UDM cosmo- 
logical model and recent observational data. Here, we 
use Cl m G [0.1, 1] and 6 X G [-1,0] in steps of 0.01. 

Figure 1 (thin dashed lines) shows the la, 2a and 
3a confidence levels in the (f2 m ,|0i|) plane when us- 
ing BAOs. Obviously, the 9\ parameter is not con- 
strained by this analysis and all the values in the interval 
— 1 < 6\ < are acceptable. However, the BAOs sta- 
tistical analysis puts constraints on the matter density 
parameter Q m ~ 0.25. 

Therefore, in order to put further constraints on 6\ we 
additionally utilize the SNIa data. In figure 1 (thick solid 
lines), we present the SNIa likelihood contours and we 
find that the best fit solution is Q m ~ 0.4 and 9\ ~ —0.05. 



The joint likelihood function peaks at f2 r 



0.2. 



:+0.02 
-0.01 



and 1 = -0.39™ (0 2 ~ 0.54) with xLtC^m, N) * 
194.1 (dof = 191). Note that the errors of the fitted 
parameters represent la uncertainties. In the insert plot 
of figure 1 we provide the solutions (circles-BAOs and 
triangles-SNIa) within la contours in the (f2 m , 4>o) plane, 
where 0o is the present value of the scalar field. The 
corresponding best fit value of the scalar field (see eqs. [15] 
and[T5l! is (j) ~ 0.42 or 0.084 in Planck units (G = c = 1), 




0.5 1 1.5 



FIG. 4: Left Panel:The equation of state parameter as a func- 
tion of the scale factor of the Universe. The lines corresponds 
to UDM (solid), CPL (long dashed) and VCG (dot dashed). 
In the insert panel we present the time evolution of T$/V(4>). 
Note, that is the kinetic energy of the scalar and V(4>) 
is the potential. Right Panel: The functional form of the 
equation of state parameter for various UDM models. The 
upper line corresponds to (Q m ,6i) = (0.73,-1) while the 
bottom line corresponds to (£l m ,9i) = (0.19,-0.1). Note 
that the solid thick line corresponds to the best-fit parame- 
ters (Sim, 61) = (0.25, -0.39). We find that initially all the 
UDM models [fi m £ (0.19,0.73), 9i G (-1,-0.1)] start from 
w — ► +1 and they reach w ~ 0(— 1) close to the present 
time. 



while the frequencies are u>i ~ 0.067Gyr~ 1 and u>2 — 
0.095Gyr~ 1 . It is interesting to mention that although 
the frequency (uj\ ~ 0.9H ) of the hyperbolic oscillator in 
the x\ axis is somewhat less than the present expansion 
rate of the universe, the u>2 is equal to the value predicted 
by the A cosmology (see section 4.1). 

Knowing now the parameter space (u)\, 9%, u)2, 02) we 
investigate, in more detail, the correspondence of the 
UDM model with the different dark energy models (see 
sections 4.1, 4.2 and 4.3) in order to show the extent to 
which they compare. Our analysis provides an evolution 
of the UDM scale factor seen in the upper panel of fig- 
ure 2 as the solid line, which closely resembles, especially 
at late times (0.6 < a < 1.5), the corresponding scale 
factor of the A (short dashed), VCG (dot dashed), and 
CPL (long dashed) . Note that the UDM deceleration pa- 
rameter at the present time is qo — —0.62. However, for 
a > 1.5, the CPL and the VCG scale factors evolve more 
rapidly than the other two models (UDM and A cosmol- 
ogy). Also it is clear that an inflection point [d(tj) = 0] 
is present in the evolution of the UDM scale factor. The 
UDM inflection point is located at tj ~ 0.46io which cor- 
responds to aj ~ 0.61 and is somewhat different than 
the value predicted from the usual A cosmology (see sec- 
tion 4.1). Before the inflection point, the UDM appears 
to be more decelerated from the other three dark en- 
ergy cosmological models due to the fact that the second 
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term [oc sinh 2 (tJit + 6\)} in ea. (fT8|) plays an important 
role. From figure 2 it becomes clear that the UDM model 
reaches a maximum deviation from the other three dark 
energy models prior to a ~ 0.15 (z ~ 5.5). In order to 
investigate whether the expansion of the observed uni- 
verse follows such a possibility we need a visible distance 
indicator (better observations) at redshifts z > 2. 

The evolution of the scalar field is presented in the 
bottom panel of figure 2, while in the insert figure we 
plot the scalar field dependence of the potential energy 
normalized to unity at the present time. As we have 
stated in section 3.1 there is one minimum at <p = 
that corresponds to t m — —0\/uj\ ~ 0At o . To conclude, 
we plot in figure 3 the relative deviations of the distance 
modulus, A(m— M), of the dark energy models used here 
from the traditional A cosmology. Notice that the open 
points represent the following deviation: (m — M)sNia — 
(to - M) A - Within the SNIa redshift range 0.016 < z < 
1.775 (0.360 < a < 0.984) the VCG distance modulus is 
close to the A one. The largest deviations of the distance 
moduli occur at redshifts around 0.5-1 for the UDM and 
1.1-1.5 for the CPL model respectively. 



4.5. The equation of state parameter 

We would like to end this section with a discussion on 
the dark energy equation of state. As we have stated al- 
ready in the introduction, there is a possibility for the 
equation of state parameter to be a function of time 
rather than a constant ratio between the pressure and 
the energy density. Within the framework of the scalar 
field cosmology the equation of state parameter is derived 
from the field model and in general it is a complicated 
function of time, even when the potential is written as a 
simple function of the scalar field. In our case we have 



w(t) 

or else (see eqsUHl [TTJ) 



P# = <j> 2 - 2V{<j>) 
P4> ~ <p 2 + 2V{4>) 



w(i) 



^ 2 (i) 



(36) 



(37) 



Note that A— models can be described by scalar models 
with w strictly equal to -1. Using our best fit parameters 
we present in left panel of figure 4, the equation of state 
parameter as a function of the scale factor for the differ- 
ent dark energy models. The UDM model (solid line) is 
the only case that provides positive values for the equa- 
tion of state parameter at early epochs. We have checked 
the UDM scenario against the cosmic coincidence prob- 
lem (why the matter energy density and the dark energy 
density are of the same order at the present epoch) by uti- 
lizing the basic tests proposed by [41|. These are: (a) at 
early enough times the equation of state parameter tends 
to its maximum value, w — > +1, which means that the 
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FIG. 5: Upper Panel: The evolution of the growth fac- 
tor for different dark energy models. The lines correspond 
to UDM (solid), VCG (dot dashed), A (short dashed) and 
CPL (long dashed) models. Bottom Panel: The deviation 
(1 — Dde / Da)% of the growth factor for various dark energy 
models with respect to the A solution. The points represent 
the comparison: (a) UDM- A (open triangles), (b) VCG- A 
(solid squares) and (c) CPL- A (open circles). 



dark energy density initially takes large values. So as 
long as the scalar field rolls down the potential energy, 
V{4>), decreases rapidly and the kinetic energy — <j) 2 /2 
takes a large value, (b) then <f> continues to roll down, the 
dark energy density decreases and the equation of state 
parameter remains close to unity for a quite long period 
of time (a < 0.2) and (c) for 0.2 < a < 0.95 the equa- 
tion of state parameter is a decreasing function of time 
and it becomes negative at a > 0.56. Before that epoch, 
the potential energy of the scalar field remains less than 
the kinetic energy (see the insert plot in the left panel 
of figure 4) and the equation of state parameter (or the 
scalar field) resembles background matter. In a special 
case where w = [or ~ V{4>)} the equation of state 
behaves exactly like that of pressure-less matter. For 
w = —1/3 we reach the same expansion as in an open 
universe, because the dark energy density evolves as a~ 2 
and has no effect on a. In fact, we verify that prior to the 
inflection point w(tj) ~ —0.334, which means that after 
ti the accelerating expansion of the Universe starts. Fi- 
nally, w ~ — 1 close to the present epoch a ~ 1 and the 
scalar field is effectively frozen (the same situation seems 
to hold also in the limit a ^S> 1). This is to be expected 
because at this period the scalar field varies slowly with 
time (see the insert panel of figure 4), so that <C V{4>) 
and the dark energy fluid asymptotically reaching the 
de-Sitter regime (cosmological constant). 
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In order to conclude this discussion, it is interesting to 
point out that we also investigate the sensitivity of the 
above results to the matter density parameter. As an ex- 
ample, in the right panel of figure 4 we present the evo- 
lution of the equation of state parameter for (f2 m ,#i) = 
(0.73,-1) [upper line] and (O m , 6>i) = (0.19,-0.1) [bot- 
tom line]. We confirm that in the range fl m S (0.19, 0.73) 
and &\ G ( — 1, —0.1) the general behavior (described be- 
fore) of the functional form of the equation of state pa- 
rameter is an intermediate case between the above lines 
for a < 1 and thus it depends weakly on the values of 
the parameter space (f2 m ,#i). Therefore, our main cos- 
mological results for the UDM scenario persist for all 
physical values of O m and it strongly indicates that the 
UDM model over-passes the cosmic coincidence problem. 



5. EVOLUTION OF MATTER 
PERTURBATIONS 

In this section we attempt to study the dynamics at 
small scales by generalizing the basic linear and non- 
linear equations which govern the behavior of the mat- 
ter perturbations within the framework of a UDM flat 
cosmology. Also we compare our predictions with those 
found for the dark energy models used in this work (see 
sections 4.1, 4.2 and 4.3). This can help us to understand 
better the theoretical expectations of the UDM model as 
well as the variants from the other dark energy models. 



5.1. The Evolution of the linear growth factor 




0.3 0.4 0.6 0.8 



t/t, 

FIG. 6: Evolution of radius of a collapsing overdense region. 
The solid and the dashed line corresponds to the UDM (S7 m = 
0.25) and A cosmology (Q m = 0.26) respectively. 

Dde(a) for the current dark energy models with respect 
to the A solution D^(a). Assuming now that clusters 
have formed prior to the epoch of Zf ~ 1.4 (af ~ 0.42), 
in which the most distant cluster has been found [69(, 
the UDM scenario (open triangles) deviates from the A 
solution by 4.2% while the CPL (open circles) and VCG 
(solid squares) deviates by —1.5% and 5.1% respectively. 
Also at the A-inflection point (af ~ 0.56) we find the 
following results: (i) UDM-A 3.3%, (ii) CPL-A -0.4% 
and (iii) VCG-A 3.2%. To conclude this discussion it is 
obvious that for a > 0.7 the UDM growth factor tends 
to the A solution (the same situation holds for the CPL 
model but with a > 0.55). 



The evolution equation of the growth factor for models 
where the dark energy fluid has a vanishing anisotropic 
stress and the matter fluid is not coupled to other matter 
species is given by ([11], [HI, [5^]): 



d 2 D 
dN 2 



l_dff\ dD_3 

H AN ) dN 2 W V ' 



where N = lna and Sl m (a) = £l m a~ 3 Hq / H 2 (a) . Use- 
ful expressions of the grow th factor can be found for the 
ACDM cosmology i n [52J| . for the quintessence scenario 
(w = const) in [55[ , |56||. |57| , [H[ , for dark energy models 
with a time varying equation of state in [59j and for the 
scalar tensor models in [60] . In the upper panel of figure 
5 we present the growth factor evolution which is derived 
by solving numerically eq. (|38|) , for the four dark energy 
models (including the UDM). Note that the growth fac- 
tors are normalized to unity at the present time. The 
behavior of the UDM growth factor (solid line) has the 
expected form, i.e. it is an increasing function of the 
scale factor. Also we find that the growth factor in the 
UDM model is almost an intermediate case between the 
VCG (dot dashed line) and CPL (long dashed line) mod- 
els respectively. In the bottom panel of figure 5 we show 
the deviation, (1 — Dde/D\)%, of the growth factors 



5.2. The spherical collapse model 

The so called spherical collapse model, which has a 
long history in cosmology, is a simple but still a fun- 
damental tool for understanding how a small spheri- 
cal patch [with radius R(t)] of homogeneous overden- 
sity forms a bound system via gravitation instability [6l| . 
From now on, we will call at the scale factor of the uni- 
verse where the overdensity reaches its maximum expan- 
sion (R — 0) and Of the scale factor in which the sphere 
virializcs, while R t and Rf the corresponding radii of the 
spherical overdensity. Note that in the spherical region, 
p mc oc R~ 3 is the matter density while p^^ will denote 
the corresponding density of the dark energy. In order to 
address the issue of the dark energy in the gravitationally 
bound systems (clusters of galaxies) we can consider the 
following assumptions: (i) clustered dark energy consid- 
ering that that the whole system virializes (matter and 
dark energy), (ii) the dark energy remains clustered but 
now only the matter virializes and (iii) the dark energy 
remains homogeneous and only the matter virializes (for 
more details see [HJ, [63| [64| and (65|). Note, that in 
this work we are using the third possibility. 

Here we review only some basic concepts of the prob- 



9 



lem based on the assumption that the dark energy com- 
ponent under a scale of galaxy clusters can be treated 
as being homogeneous: p<p c (t) = p<f,(t), </> c (i) = and 
w c (t) — w(t). In general the evolution of the spherical 
perturbations as the latter decouple from the background 
expansion is given by the Raychaudhuri equation: 

3R = -4nGR[p mc +p< j>e (l+w c )] here 4ttG=1/2. (39) 

Now within the cluster region the evolution of the dark 
energy component is written as (see [62]) 



R 

p<t>c + 3 ^( 1 + W0 c )p0 c =r 



(40) 



while if we consider a scalar field the above equation be- 
comes 



3^c + t/'( 



where 




(41) 



(42) 



Figure 6 presents examples of R(t) obtained for the UDM 
(solid line) and for the concordance A model (dashed 
line). The time needed for a spherical shell to re-collapse 
is twice the turn-around time, tf ~ 2t t . 

On the other hand, utilizing both the virial theorem 
and the energy conservation we reach to the following 
condition: 



\r-^(u g + u iK ) + u g + u^ 



[Ug + U* B 



(43) 

where Ug = —3GM 2 /5R is the potential energy and 
U<i> c = -4nGM(l + 3w c )p0 c i? 2 /5 is the potential en- 
ergy associated with the dark energy for the spherical 
overdensity (see [62[ and [H[; in our case 4itG = 1/2). 
Using the above formulation we can obtain a cubic equa- 
tion that relates the ratio between the virial i?f and the 
turn-around outer radius i?t the so called collapse factor 
(A = Rf/Rt). Notice that eq. (|43|) is valid when the ratio 
of the system's dark energy to the matter's densities at 
the time of the turn-around takes relatively small values 
(64j. Of course in the case of w c = —1 the above ex- 
pressions get the usual form for A cosmology ([13], (66|) 
while for an Einstein-de Sitter model (tt m = 1) we have 
A = 1/2. Finally solving numerically eq. (|4"3")) [it can be 
done also analytically] we calculate the collapse factor. In 
particular, figure 7 shows the behavior of the collapse fac- 
tor for the current cosmological models starting from the 
UDM (solid line), A (short dashed), VCG (dot dashed), 
and CPL (long dashed) . We find that the collapse factor 
lies in the range 0.43 < A < 0.50 in agreement with previ- 
ous studies ( [H3, [U [U, [13, [13, HI, [H). Prior to 
the cluster formation epoch (zf ~ 1.4) the UDM scenario, 
appears to produce more bound systems with respect to 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 


1 1 








MUM 

A-Cosmology 

7 CPL " 

! VCG : 

" i 1 i i i 1 i i i 1 i i i 1 , , i 1 i i i 1 i " 


0.4 0.6 0.8 1 


1.2 


1.4 



FIG. 7: The collapse factor versus the redshift of virialization 
for various dark energy models. 



the other dark energy models. Indeed, we find the fol- 
lowing values: Xudm — 0.44, Aa ~ 0.49, Xcpl — 0.48 
and Xvcg — 0.50. Also it becomes clear that the UDM 
collapse factor decreases slowly with the redshift of virial- 
ization Zf, due to its positive equation of state parameter. 
This is also incorporated by the fact that at early epochs 
the cosmic expansion of the UDM model is much more 
decelerated than in the other three dark energy models. 
The latter result is in agreement with those obtained by 
[62| . They found a similar behavior for the collapse fac- 
tor by considering several potentials with an exponential 
phase. 



6. CONCLUSIONS 

In this work we investigate analytically and numeri- 
cally the large and small scale dynamics of the scalar 
field FLRW flat cosmologies in the framework of the so 
called Unified Dark Matter scenario. In particular using 
a Hamiltonian formulation we find that the time evolu- 
tion of the basic cosmological functions are described in 
terms of hyperbolic functions. This theoretical approach 
yields analytical solutions which can accommodate a late 
time accelerated expansion, equivalent to either the dark 
energy or the standard A models. Furthermore, based on 
a joint likelihood analysis using the SNIa data and the 
Baryonic Acoustic Oscillations, we put tight constraints 
on the main cosmological parameters of the UDM cos- 
mological model. In particular we find Q m ~ 0.25 and 
the scalar field at the present time is cf>o ~ 0.42 or 0.084 
(in Planck units). Also, we compare the UDM scenario 
with various dark energy models namely A cosmology, 
parametric dark energy model and variable Chaplygin 
gas. We find that the cosmological behavior of the UDM 
scalar field model is in a good agreement, especially after 
the inflection point, with those predicted by the above 
dark energy models although there are some differences 
especially at early epochs. In particular, we reveal that 
the UDM scalar field cosmology has three important dif- 
ferences over the other three dark energy models consid- 
ered: 

• It can pick up positive values of the equation of 
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state parameter at large redshifts (z > 0.8). Also, 
it behaves relatively well with respect to the cosmic 
coincidence problem. 

• At early enough epochs (a ~ 0.15 or z ~ 5.5) 
the cosmic expansion in the UDM model is much 
more decelerated than in the other three dark en- 
ergy models. In order to investigate whether the 
expansion of the observed universe has the above 
property, we need a visible distance indicator (bet- 
ter observations) at high redshifts (2 < z < 6). 

• Close to the cluster formation epoch, its collapse 
factor Xjjdm is less than 12% of the corresponding 



factor of the other three dark energy models. This 
feature points to the direction that perhaps the A 
parameter can be used as a cosmological tool. 
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